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ABSTRACT 


This thesis addresses the problem of estimating the distance which 
must be traveled to visit all the targets in a search area. The search 
area is assumed to be a compact area drawn from an infinite Poisson 
field. Targets are visited in "nearest uninvestigated neighbor" order. 
A simplified analytical model is developed for the unbounded 
investigation. This model is used to provide heuristic justification 


for simulation results obtained for the bounded investigation. 
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I. INTRODUCTION 


A. GENERAL DESCRIPTION OF THE PROBLEM 

This thesis addressed the problem of estimating the distance 
which must be traveled to visit targets which are randomly distributed 
over a large area. The problem is interesting because many practical 
search operations require the investigation of targets after detection. 
It is often necessary to estimate the time which will be required for 
the entire operation before the targets are even detected. 

Classical search theory addresses cases where it can be assumed 
that there are few targets in the area, and that all are of equal 
interest when detected. Later classification may rule out some targets 
as false targets. This is not usually obvious to the searcher. Any 
detection will stop the search. Subsequent classification as a false 
alarm may result in its restart. 

A variation on this case occurs when there are many false targets, 
but they may be classified immediately upon detection. If the target 
may be immediately classified as false, then the search can proceed 
unimpeded. False targets can be treated as if they had never been 
detected. If all false targets are immediately classified as false, 
then this case becomes equivalent to the classical case of sparse 
targets. The search continues until the true target is found or the 
area is exhausted. When the first target is detected which cannot be 
dismissed as false, the search is stopped and that target is pursued. 

There is a third case which shall be addressed. It is inter- 


mediate between the first two. It can be defined as follows: targets 





can be easily classified, but only at very close range. The searcher 
can detect targets at long range, but, for classification, must take 
a "close look" at which time he knows immediately whether or not he 
has found the object of his search. If he has, he stops, if he hasn't, 
he keeps searching. If the. targets are dense enough, or the searcher's 
detection range great enough, the searcher may always have uninvestigated 
targets in sight. He may not really "search" at all in the sense of 
classical theory; he may instead travel directly from target to target 
without any regular "search pattern." The problem is now one of 
"investigation theory" as first defined by Nuendorfer [1]. The 
investigation will continue until the desired object has been found, 
or until all of the targets have been visited. In either event, the 
investigator's track will have gone to the location of each prior 
false alarm. 

This is the crux of the problem. If the positions of the targets 
are unknown until detection, then the length of the entire path will 
be uncertain. If the target positions are random variables, then the 
length of the path required to investigate the targets will also be 


a random variable. 


Be EXAPIPLE 

The case described above is typical of many Coast Guard operations. 
The Coast Guard is often required to search for a particular boat. 
Usually the boat can be positively identified only by name or registra- 
tion number. Small boats which have been reported overdue, and boats 
suspected of law enforcement violations are two classes which must be 
searched for in this manner. There may be many boats which fit a 


general description, and the searcher may have to check out all of them. 





A simple scenario can be constructed which illustrates the problem. 
Consider a Coast Guard Cutter with embarked helicopter assigned to 
patrol a ten thousand square mile area for fishing violations. The 
helicopter must be fueled tonight for tomorrow's first light patrol. 
The aircraft commander must decide how much fuel to load. Too little 
fuel means an incomplete patrol. Too much fuel means that something 
else must be left behind. Experience has shown that rarely are more 
than twelve or fewer than eight ships seen. Ten is the average number. 
The aircraft commander has no way of predicting the future position 
of any of the fishing vessels, but he must still make an important 
operational decision, based on the length of the path which he expects 
to fly. The result developed in this paper would be directly 


applicable in this scenario. 





II. ANALYSIS 


A. ANALYTICAL FORMULATION OF THE PROBLEM 

Two basic questions must be addressed first off: 

1. What distribution is to be assumed for the positions of the 
targets? 

2. In what order are the targets to be visited? 

In this thesis it is assumed that the area considered is a portion 
of an infinite Poisson field with targets distributed with spatial 
Poisson density p. This means that the number of targets in any 
finite area, A, will be a random variable, with a Poisson distribution 
and mean pA. P.A.W. Lewis and G.S. Shedler [2] have shown that the 
targets occurring in an area as the result of this process will be 
uniformly distributed within the area. This is to say that if the 
entire area is divided up into many small boxes, each large enough to 
hold only one target, then each box has an equal chance of actually 
holding one target. This distributional assumption is attractive 
because it assumes the minimum amount of information about the target's 
location. In this sense it is the “worst case"; more information can 
only improve the situation. In actual Coast Guard practice, there are 
many searches conducted when uniformity is the best assumption which 
can be made; the actual distribution of targets, including false targets, 
cannot be distinguished from what would have been expected from a 
uniform distribution. 

A related question is that of target motion. In this thesis it 
will be assumed that the targets move so much more slowly than the 


investigator that they may be considered stationary. 





The question of the order in which targets are to be visited is 
more complicated. If the targets are stationary, there is one unique 
optimum order resulting in the minimum total path length. The problem 
of determining this optimum ordering is called the "travelling 
salesman problem." The travelling salesman solution, while ideal, 
is not practical for large problems. There is no efficient method of 
solution. Every current method is fundamentally equivalent to trying 
all possible orders [3]. There are some shortcuts for smal] numbers 
of points; the fact that points on the perimeter will be taken in the 
relative order that they are found around the perimeter is an example. 
But these shortcuts do not work for large numbers of points with 
complex geometry. Determining the optimum path through twenty points 
will require a tremendous amount of calculation. 

This thesis employs the myopic decision rule in which the investi- 
gator will always go to the nearest point which has not yet been 
visited. This decision rule is justified on two grounds: 

1. It is easily employed in actual searches; the searcher need 
only remember what targets have been investigated, and then select the 
closest of the remainder. 

2. Balut [4] has shown simulation results indicating that the 
path resulting form the myopic decision rule 1s only about 15% longer 
than the optimal path for investigations where the investigator is 
much faster than the targets. This procedure is often used in actual 
searches. 

Finally, the assumption is made that the investigator's detection 
range is much greater than the inter-target spacing. This means that 
the investigator will always know the location of the nearest uninvesti- 


gated target. 





B. MODEL FOR AN UNBOUNDED MYOPIC INVESTIGATION 

Because the area considered is a portion of an infinite Poisson 
field, it is profitable to first consider the properties of an investi- 
gation conducted in an infinite field. The investigator will start 
at one of the points in the field and then visit its closest neighbor, 
then the closest neighbor remaining, and so on. The investigator will 
always travel to the nearest uninvestigated target (n.u.t.). Real 
searches never detect an infinite number of targets, but this is not 
a bad approximation to the case when the area is very large and there 
are many more targets detected than can be investigated in the course 
of a mission. 

Consider a Poisson field with spatial density p. Call the set of 
points I = {Py> Py> Pos Pas veers P,? an investigation starting at Po: Iie 
Py > Pop j.e. the investigator returns to the point where he started, 
then the investigation is a tour. The investigation is conducted 
according to the following rule: the investigator will always visit 


the closest point to which he has not yet been. Let r be the distance 


from Psy to p.. Let dy = fore Let ba = 


‘ rss when Pn = Po- 


irda 
| 

tois 
paar 


by is the cumulative distance convered in visitina the first k noints; 
to is the length of a tour visiting all the points and returning to 
its starting point. 
The expected length of the first radius, ry> may be calculated as 
follows. The probability of a target being found in a small area, dA, jis 
P(dA) = pdA. 
Thus the probability that a target will not be found is 
P(dA) = 1 - pdA. 


1] 





Let F(A) be the probability that there is no point in some large area A. 
The probability that there is no point in (A + dA) is 
F(A) + P(dA) = F(A)(1 - pdA). 
But this is also F(A + dA) 
'. F(A)(1 = pdA) = F(A + dA) 
thus 


Fins = - F(A) _ “pF (A). 


As dA + 0, then this becomes 


dF(A) _ = 

“da PFA) 
so that 

F(A) = ce Ph 


where c is the arbitrary constant of integration. 
feo) = 1, then c = 1 and 
F(A) = ae 

The area considered so far has been perfectly general, of arbitrary 
shape. 

Now consider A to be the area circumscribed by a circle of radius r, 
A = ar’. Let A, be the smallest such circle tangent to an adjacent 
point. Let r, be its radius. The tail distribution of A, Fa, (t) is 
defined as follows: 


Fy (t) = P(A, > ¢) 


1 ] 


Let t = We. where jt is the area circumscribed by a circle of radius x. 
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Then 
Fa (t) = P(A, > t) 


] 
= P(nrs > eel 
= P(r,“ > x5) 
= P(r, Sas X»0y es 0) 
Define 
8.) = P(r, > x) 
= fF C3) 
7 
But 
E Gee ePt 
“7 2 
-Pp1rXx 
2 
Cae oa 
a 


The expected value of ry E(r,); is defined as 


“ 
wef jedi u du 


Integrating by parts yields 


E(r, ) = |, ss an = A exp [-pru“] du 


where G. Coe P(r, > x) 


Spatial Poisson density. 


ao) 
it 


ie 





Thus the expected radius of the smallest circle containing an adjacent 


Beant 1s 


Figure I[I[-] 


The expected values of successive jumps are much more difficult 


to determine. The distribution in the fresh (unsearched) area is still 


but the new area is no longer a circle. Let Xy be the actual distance 
between Po and its nearest neighbor, Py - Then a circle of radius xy: 
centered at Po: is Known to contain no other points besides Po: The 

new point, Py: lies on to this circle. Area included in this circle 
“doesn't count” in the search for new points and must be ignored in 
calculating the distribution function for the next interpoint distance, 
Geetx). Let H. (x) be the function whose value is the new (uninvestigated) 


a 
2 
area lying within a circle of radius x, centered at Pie 











2 
Area known to be 
free of points 
Figure II-2 
Here aes ers, thus 
r c if r ie c 
ae ey) = (154 - r)° ci : 7 - 1) cos | = - = a 
or ] ] Ar 


Se? 


im - a at 
FO a ae ; : 
. c © J 7 a 





Thus 


Hy) (x) = nxe and G (x) =e 


It can easily be seen that 


Hy (x) = x(x° - x,°) 


2 Z 
2 2 X -] x x x 
Ske = Xs = (2x. ( Joos eo - = V1- } 
1 ] ae 2x, 2x, ay 2 








ix < 2x, ) where x, is the actual observed distance between Py and p, 


Let G(x) = P(r, exe 


Then, as before 
-pH, (x) 


and 





2x, 9 
| exp{-p{7(x - x°) - [2x, ( : i 1 cos” ae = 
0 


2X, 
/ aye 5 ]}}dx 


4X, 


| exp { -pr(x¢ - x1") }dx 


2x, 





This expression is difficult to evaluate analytically, and the expected 
values of ra, Ty, ... are even more difficult to obtain. This is because 
the expected value for re depends not only upon ry and Po see Te lys 

but also upon the angles between successive radii. An investigation 


which doubled back on itself, Figure II-3(a), would have a different next 


expected radius than one which did not. 
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ae) 
~_— 
ar 
— 


Figure [[-3 


While the total sanitized area is the same in both cases, it is 
not distributed according to a common function of radius. 

Because the distributions and moments are so difficult, approxima- 
tions would be useful. An argument can be advanced to support the 
conjecture that E(r.), i wieeeee eee are asymptouically |inear as 
1 +o, Because new points will be located in an uninvestigated area, 
the process will tend to creep away from sanitized areas and toward 
unsanitized areas. This will usually put the investigator at the 
boundary of a fresh area, with a constantly expanding sanitized area 
benind him. Asymptotic linearity may be shown easily in a restricted 
model based on this behavior. Consider a completely deterministic 
process where points appear on a straight line. The spacing between 


the points will be determined as follows: 


ee 





lz ry will be arbitrary, 


a r. will be such that 





meeeall in i> I. 

This is equivalent to circumscribing a circle of radius ry about the 
first point, and then adding other points so that the new area added 
to the process at each step is equal to the initial area enclosed in 
the first circle. The process can be visualized as a sequence of 
circles proceeding down a line. Each new circle is centered on the 
intersection of its predecessor and the centerline. Each circle 1s 
sized so that its area, less the area of its predecessors, iS equal to 


ry’. Refer to Figure I[I-4. 


Figure I[-4 


Numerical calculation for all initial values selected between 


1 90 


i walla 10° and r, = 1.0 x 10°" have shown that the sequence 


Tyo Pos Vas sees V . 1S asymptotic to 1.2814 ry. 


n? ee 
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While one would never expect to see an actual realization develop 
along a perfectly straight line, an extensive computer simulation by 
Laidlaw [5] does show that many investigations do develop in a way that 
could be characterized as random deviations on either side of a mean 
line of progress. It seems reasonable that, if this were so, and the 
investigation did proceed as variations about a mean line, then the 
actual limiting value might be r, > cr), where c > 1.2814. Laidlaw's 
Simulation yielded c) = 1.405 based ona myopic investigation conducted 
in a finite plane containing 10,000 points, with the investigation 
stopping when the boundary of the plane was closer than the nearest 
target, and with the sweep width much greater than the expected target 
spacing. 

The assumptions for this result were: 

1. The points visited are located in an infinite Poisson field. 

2. Detection range is much greater than target spacing. 

3. Myopic search is employed. 

Under these assumptions the length of the total path appears to 
be asymptotically convergent to 


2 =ner 
n 1 


n(1.405) 
2vD 


The model is asymptotically linear in n and is called the asymptotic 


linear model of the unbounded investigation. 





C. APPLICATION TO THE BOUNDED INVESTIGATION 

The linear asymptotic model cannot be directly applied to the 
bounded investigation. The linear model requires that the spatial 
Poisson paramenter, Pp, remain constant in the unsearched area. This 
will not be true ina finite area, because the finite area will 
contain only a finite number of points. 

Nevertheless, the linear asymptotic model does suggest an indirect 
approach. Consider a bounded investigation to consist of two parts: 

1. An investigation of the points near the center of the area 
where the unbounded linear model may be applied. 

2. An investigation of the remaining points near the boundary. 
It seems reasonable to suppose that the bounded investigation might 
first behave like the unbounded case, when there are a large number 
of points from which to draw. After a while the points will have 
been "thinned out" considerably, and the investigator's path will take 
him back to "mop up" the remainder. 

A simulation was constructed to test this conjecture. A number 
of points were generated in a square. An investigation was started 
at a point near the center. The path took the investigator to all 
the points in the area, and then back to his initial point. Such a 
closed path is called a tour. Plots of 2) /KE(Y, ) - vs - k for 
meinen, 1 = 10, 30), 60 shown in Figure I1I-5. A cumulative 
Dot Of t/nE(r,) - vs - n is shown in Figure III-6. 

Two things can be seen in these graphs: 

1. The “plateau" in the center of the curves in Figure III-5 becomes 
more pronounced as more points are added to the area. This 1s a 


reasonable result. 
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K, THE NUMBER OF POINTS VISITED SO FAR. K = 1,2,...4 


HATIO OF CUMULATIVE PATH LENGTH, Le 180, 


Figure [I-5 


2] 


(KxECR ,1)=VS=K,K=1,2,...5M 





RATIO: (T,,)#(WxECR,]) 
2.16 __ 


NV, THE TOTAL NUMBER OF POINTS IN THE PLANE 


RATIO OF TOTAL TOUR LENGTH, Ty, TO (NxECR,1)=VS= W 


Figure [[-6 
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2. The final cumulative length is approximately 


t = 1.9508) = 1.95nE(r,), hee 


Plots of repeated simulation trials suggest that t, 1s normally 


distributed with mean 


- _ nA 
£ = 1.95 3A 
s * = 9.06 TA 
t Z. 


Stemleaf plots for the cases n = 20, 40, and 60 are shown in 
pugure III-/. 

This is the expression for t when n is known. Generally n will 
not be known, but the assumption that the search area is taken from 
a Poisson field will Baan the calculation of 2105) from p and A, 


The expected number of points occurring in an area A is E(n) = pA. 





Thus 
E(t |") = 1.95uE(r, ) 
= 508 tee 
2v p 
and 


y 1.95 k LL F£ (k) 
k=0 DA 


7 
-_— 
ct 
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Similarly, 


3/2 
Sc, (0y {05 ae 
t. q 


III. EXAMPLE 


The procedure developed may be applied to the example given in 
seetren II. 

Assume that the search area is a square area, 100 miles on a side, 
with a spatial Poisson parameter of 


ye i 6 Wore 


The expected length of a tour will be given by E(t, ) =| 28)s tc 
thus 
E(t) = 323 miles. 
To check this result, eleven points were selected at random, with 
one being the searcher. The positions, measured from the southwest 


Gomer, were: 


Index Xx y 
] Sl Sir. 
2 60.5 00.4 
5 90.5 ges 

49.8 em 
5 18.4 Sieh 2 
6 65.3 WES, 
7 40.8 Zee 
8 51.6 Ly died 
9 69.7 99.3 
10 58.0 74.0 
0 18.0 19.8. (searcher) 
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These points are plotted in FigureIII-I. It is assumed that the 
searcher's detection range is greater than the greatest inter-target 
spacing. 

The myopic investigation starting at point "0" is shown by the 
memceiane in Figureiii-~2. [I = (0, 7, 4, 8, 2, 3, 6, 10, 9, 5, 1, O}. 
The actual length of the myopic investigation starting at pt 0 is 
328.66. In this case the myopic path is also the optimal path through 
all the points, and 328.66 1s the minimum length. In this case n = 11 
and 
s, 2 = .06 1650 


an 


Saas 40-62 
th 


So the actual value 1s well within one standard deviation of the 


predicted value. 
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Figure III-] 


Point Locations for Example Problem 


zi 








Figure [[I-2 


Myopic Tour Through Example Problem 
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IV. 


SUMMARY AND CONCLUSIONS 


This thesis has presented a method of estimating the length of a 


path required to visit all of the targets in a search area under three 


conditions: 


|. The targets are distributed uniformly within the area. 


2. The "nearest uninvestigated target" (myopic) decision rule is 


used. 

3. Detection radius 
Simulations results have 
tour of all points under 
a simple function of the 
is known, or by the area 


targets is not known. 


1s much greater than inter-target spacing. 
Shown that the expected distance of a complete 
these assumptions is predicted very well by 
area and number of targets, if that number 


and average target density, if the number of 
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RAT IOS<0 

OMS "CLOSIO PRINTER OFF' 


? 


DD D D 


» 


OFFLINE PRINTING OF INTERMEDIATE AND FINAL RESULTS IS CONT 


t 


'"VUMBER OF TRIALS THIS RUN '3NTRIALS 
"NUMBER OF POINTS GENERATED PER TRIAL 
"SPATIAL POISSON CONSTANT, LAMBDA 
GX LENGTH 2x LAMBDA*0.5 

"EXPECTED LENGTH OF FIRST JUMP ';&XLENGTH 
? 


SV EOlVe o. 
' ; LAMBDA*NPOLINTS*+10*6 


? 


CUMS+0 
INDEX*NTRIALS 


Sree tiesPOlNTS 


B<( WPOINTS?1000) ,C1.5](NPOINTS?1000) 
COUNT+1+0B 


STEP2:A~, BLCOUNT ; ] 


>STEP3x1(502] 500=4[1])x502]500=A[2] 
>STEP2* COUNT +COUNT =1 
+STEP 


STEP3:B+ 10 +(COUNT=1)08 


RATIOS*RATIOS,A NUT2 B 

CUMS+CUMS+0U TPUT ,C2](OUTPUTL 36)*2) ,01.5](OUTPUTE 37]*2) 
+STEPX\INDEX*INDEX=1 

CUMS*CUMS*NTRI ALS 
CUMSL:8]<(NTRIALS+NTRIALS=1)x(CUMSL 381]=CUMSL 36)*2) 
CUMSL:9]<(NTRIALS+NTRIALS =1) x(CUMSC39]=CUMSL ;71*2) 
CH+NUT2CH,',L VAR,R VAR' 

HATRIALSFMT FMT CUMS 

+NOPFLINEX\OPRINT=0 

OPRINT A 

'CLOSIO PRINTER ON' 

RATIOS+1+ RATIOS 
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ff ) 

2 ] 

ea 

C4] 

C5] 

[6] 

7) 

[8] 

C9] 

[10] 
hai) 
in } 
[13] 
C14] 
C15] 
[16] 
7a 
[138] 
C19] 
[20] 
C21] 
22 } 
(23) 
[24] 
[25] 
C26] 
27 | 
[28] 
[29] 
[30] 
C31] 
aval 


C33] 
C34) 
C35] 
[36] 


V 


V 


VWUT2COIV 
PAT4HLTA*STARTSTOP NUT2 POINTS; DIST ;DIST2;UNVISITED;NEXT ;NOW 3NUM; 
NUM213START ;STOP 
a TaIS FUNCTION FINDS THE MYOPIC INVESTIGATION THROUGH THE 
POINT SET 'POINTS' STARTING AT 2+STARTSTOP AND ENDING AT 
“Q+STARTSTOP. OUTPUT IS THROUGH THE ARRAY ‘OUTPUT! AND 
CONSISTS OF THE LIST OF POINT INDICES (FROM 'POINTS'), POINT 
COORDINATES, JUMP LENGTHS, CUMULATIVE LENGTHS, AND RATIOS 
OF CUMULATIVE LENGTH TO (NxEXLENGTH), WHERE 'BXLENGTH' IS THE 
a EXPECTED LENGTH OF THE FIRST JUMP. 
NUM+SIZE 
VUTJUMPS+*0 
VUTLENGTHS+*0 
LENGTH<0 
NUTPATH<NOW+1 
LIST2+(NUM+1) +LIST+1N+2+1+p POINTS 
START+2+STARTSTOP 
STOP*+ 2+STARTSTOP 
DIST+DISTANCE START ,{1] POINTS,{1] STOP 
DIST+DIST+ZAPx(1N)°.=( 14) 
UNVISITED*0,((N=2)01),0 
STEP :DIST2<UNVISITED/DIST 
NEXT+1+A, DIST2(NOW ; ] 
LENGTH<LENGTH+DUMMY+«, DIST2( WOW 3 VEXT] 
VUTIUMPS*VUTIUMPS , DUMMY 
WUT LENGTYS*NUTLENGTHS , LENGTH 
NOW+14 (NEXT =1)OUUVISITED/ LIST 
UNVISITEDLNOW J+0 
NUTPATY*NUTPATH , NOW 
>+STEPX YOZNUMANU = 1 
NUTPOINTS<START ,(1] POINTS(1+WUTPATH=13;1,01] STOP 
VUTPATA+( NUTPATH=1) ,0 
NUTLENGTHS<NUTLENGTHS , LENGTH+PATHLTH+, DISTLNOW 3N J 
NUTJUMPS+NUTJUMPS , PATHLTHY 
OUTPUT+( LIST2,1),£2] WUTPATH,(£2) NUTPOINTS,(C2] NUTJUMPS,(2] 
NUTLENGTHS ,€1.5](NUTLENGTHS+EXLENGTHx( .oNUTLENGTHS ) =1) 
PATHLTH*OUTPUT( NPOINTS+137] 
+OFFLINExQPRINT22 
OFFLINE :CH+NUT2CH 
OPRINT NUT2FPMT FMT OUTPUT 


DDD DPD D 
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VDISTANCELUIV 
V DIST*DISTANCE POINTS 
mi | A THIS FUNCTION GENERATES AN NxN DISTANCE MATRIX 
M2) DiST<(POINTSL;2)°.=POINTSL;2])+*2 
Mey) «PEST <DIST+(POTWTSL:31J°.=POINTSL31))*2 
mt) olST<DIST+0.5 
V 
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